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In this paper a new class of formations that maintain a constant 
shape as viewed from the Earth is introduced. An algorithm 
is developed to place n spacecraft in a constant shape forma- 
tion spaced equally in time using the classical orbital elements. 

To first order, the dimensions of the formation are shown to be 
simple functions of orbit eccentricity and inclination. The per- 
formance of the formation is investigated over a Keplerian orbit 
using a performance measure based on a weighted average of the 
angular separations between spacecraft in formation. Analytic 
approximations are developed that yield optimum configurations 
for different values of n. The analytic approximations are shown 
to be in excellent agreement with the exact solutions. 

INTRODUCTION 

Clusters of low-performance spacecraft flying in formation may provide enhanced 
performance over single high-performance spacecraft. This is especially true for re- 
mote sensing missions where interferometry or stereographic imaging may provide 
higher resolution data. There are many possibilities for configuring formations of 
satellites for such missions, but there are relatively few established performance mea- 
sures for evaluating the effectiveness of particular configurations. This is especially 
important as configurations vary during an orbit due to orbital dynamics, and over 
longer time scales due to perturbations. 

Significant attention has been given in the literature to many aspects of the dis- 
tributed spacecraft concept. A preliminary feasibility study of formation flying tech- 
nologies was performed by Folta et a/. 1 At the time of study it was determined that 
limitations in pointing accuracy could not meet an 80% overlap in the presence of 
errors. Further investigations were performed by Folta et al. 2 to understand how to 
relate observations from two spacecraft, to determine the effects of perturbations, 
and to confirm previous results. The formation flying concept was revisited by Folta 
et al. 3 in 1996 and the feasibility of formation flying for particular missions in light 
of new advances in technology and autonomy was investigated. The dynamics of an 
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uncontrolled formation with and without atmospheric drag was studied. A formation 
control strategy for the Earth Observing System-AMl (EOS-AMl) and Landsat-/ 
was investigated using a dynamic ballistic coefficient strategy. 

DeCou 1 developed a three-satellite station-keeping strategy describing paths to 
perform very long baseline optical interferometry. DeCou concluded that the thrust- 
ing required is low enough to allow missions of up to ten years using ion thrusters 
without refueling. Gramling et al . 5 demonstrated the possibility of using Goddard 
Space Flight Center’s (GSFC) Onboard Navigation System (ON'S) for relative nav- 
igation of the EO-l/Landsat-7 spacecraft formation. Hartman et al 6 extended the 
study of GSFC’s ON'S performance for spacecraft formation flying. Specifically Hart- 
man investigated orbit determination and prediction accuracies and their impact on 
the relative separation errors. 

The use of GPS for relative navigation of formation flying spacecraft has been 
studied by Guinn et al. 7 and for both relative position and attitude sensing by Adams 
et al. 9 and How et al. 9 Adams et al. 9 presented results from a GPS-based relative 
navigation and attitude sensing system developed in the laboratory. 

In this paper, we introduce a class of formations we call rotating formations, where 
the configuration remains relatively constant during an orbit. A detailed analysis of 
the performance of this class of formations, including an approximate analytical de- 
termination of the optimal configuration based on an angular performance measure is 
presented. The new class of formations involves an elliptical arrangement of satellites 
in slightly elliptical orbits, each inclined slightly with respect to the orbital plane 
of the reference circular orbit. Each orbit has the same semimajor axis, inclination 
with respect to the reference plane, eccentricity, and argument of periapsis. The 
position of each satellite is determined by these common elements, and by unique 
values for the right ascension of the ascending node and true anomaly at epoch. A 
general expression for these angles is developed in terms of the number of satellites 
in the formation, n, based on a equal-time distribution of the satellites. The relative 
configuration remains nearly constant during an orbit, as the satellites rotate about 
the reference circular orbit. In our development, the eccentricity e and inclination 
i are treated as design variables, and values are chosen to optimize an appropriate 
performance measure. 

One measure of the performance of such a formation is the angular separation 
between the satellites, as this defines the performance of interferometry and of stere- 
ographic imaging. A particular instrument s effectiveness in these missions can be re- 
lated to the angular separation, a tJ , between the i t h and j t ,h satellites in the formation. 
For a specific value, a m , the performance of any pair of satellites will be maximized. 
Since the angular separation varies during the orbit, the performance measure is ac- 
tually an orbit-averaged quantity depending on all the angular separations. For the 
rotating configuration developed here, the integrand of this performance measure is 
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neatly constant. especially for large n, since the time it takes the formation to reach 
a nearly ('((iiivalent conligtiration is T / n, where T is the orbital period. Tims the 
integrand can be used as a proxy performance measure, which allows the analytical 
development of the measure with all orbital elements included in the (expression. 

The full analytical expression of the performance measure is extremely compli- 
cated, and its complexity increases geometrically with n. We develop a relationship 
between the eccentricity and inclination that permits us to develop a series approx- 
imation for the performance measure. Using only terms to second order, simple 
approximations are developed for the values of e and i that maximize the formation 
performance. These approximations take the form of a function that depends only 
on n multiplied by the optimal separation angle a m . Thus the optimal e and i can 
easily be computed for the different optimal separation angles associated with dif- 
ferent instruments. The approximations exhibit excellent agreement with the exact 
solutions obtained using the orbit-averaged “exact” version of the performance mea- 
sure. Furthermore, the functional dependence of the optimal e and i on the number 
of satellites, n, is recaptured by investigating the solutions for different values of n. 

DEVELOPMENT OF THE ROTATING FORMATION 

In this section, a formation of n spacecraft is developed such that the satellites 
move along an ellipse-like shape in a relative reference frame. A configuration of n 
spacecraft on an ellipse-like path is shown in Figure 1. In general, the horizontal 
dimension is defined as the longitudinal separation and the vertical separation as the 
latitudinal separation. The dimensions can be physically represented using the angu- 
lar or spatial separations. To distinguish between the two representations, the angular 
dimensions of the ellipse are measured using the angles 6 i m and 5i at to specify the 
longitudinal angular separation and the latitudinal separation respectively. When the 
dimensions of the ellipse are discussed in units of distance, the horizontal dimension 
is termed di on and the vertical dimension is called dt at - 

In the development of the rotating formation two coordinate systems are used: an 
inertial frame Tr, and a relative reference frame JF re , defined by a reference circular 
orbit. The reference circular orbit is chosen to be equatorial with a true longitude 
at epoch of i re f. The origin of T re i is attached to a reference spacecraft that moves 
along the reference orbit. The x r6 /-axis is in the direction of the reference spacecraft’s 
position vector. The t/ re /-axis lies along the reference spacecraft’s velocity vector 
direction. Bv this definition the y re i- axis is always tangent to the reference orbit. 
The 2 > e /-axis is chosen to complete the right-handed set, and is in the orbit normal 
direction. 

Frame Tr is an inertial frame but it is not equivalent to the Earth-centered inertial 
frame T\ commonly found in the literature. The frame is fixed in inertial space with its 
origin at the center of the Earth. However, its axes can be defined by any arbitrary 
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Figure 1: Geometry of the Rotating Formation 

right-handed system. The result is that the xr-ur plane of Tr is not the Earth’s 
equatorial plane, it is the orbital plane of the reference circular orbit. The orbital 
elements used to describe the rotating formation are used according to their classical 
definition, but for generality they are measured with respect to Tr unless otherwise 
stated. 

We wish to determine the dimensions of the relative orbit in terms of the orbital 
elements of the spacecraft in formation. To simplify the analysis the discussion is 
limited to orbits that are nearly circular and result in small displacements from the 
reference orbit. In doing so, we require that the spacecraft motion occurs in nearly 
the same plane as the reference motion. Therefore, the along-track displacements in 
the relative motion are primarily due to the eccentricity of the orbits. With these 
assumptions the dependence of 5[ on and 5i at on the orbital elements can be determined. 

A view of the ellipse-like path of the rotating formation in Tr is shown in Figure 2. 
The vector r x describes the position of an arbitrary spacecraft. The true anomaly, u, 
for the spacecraft is measured from the i x axis. The quantity yi m is the along-track 
component of angular separation of the spacecraft in T re i\ he., yion is defined as an 
angle and not a distance. From inspection of Figure 2, we see that the angle 8 lat is 
due to the difference in inclination between the orbit and the reference orbit. The 
expression for 8i at in terms of the inclination is simply 

6lat = '2l (1) 

Recall that it is assumed that the out-of-plane motion does not affect the along- 
track angular separation. Therefore, an expression for S lon can be determined by 
investigating the projection of the motion in the xr — yR plane as seen in Figure 3. 
First an expression for yi on is developed. From inspection of Figure 3. y Xon can be 
written 

Vum = v-9 (2) 
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Figure 2: The Rotating Formation Seen in J-r and T re i 


where 6 is the angle between i re i and ij. An expression for 9 can be written as 


6 — M 0rel + Tjt 


(3) 

where Mo rel is the angle between i re / and ii at the initial epoch. The horizontal 
dimension, 5^, is simply twice the maximum in the along-track separation: 

&lon — 2yion\max 


(4) 

For small e we can write 



v = Mq 1 + Tjt + 2e sin (M 0l + rjt) 


(5) 

as found in Battin, 10 where is the mean anomaly of iq at the initial epoch, 

77 is the mean motion. Substituting Eqs. (3) and (5) into Eq. (2) we obtain 

and 

Vion — 2e sin (Mq x 4- r]t) + Mq x — Mq tc1 


(6) 

Recall that Mq x and Mq tc1 are constants so 



Vlonlmax 2e 


(7) 

Substituting Eq. (7) into Eq. (4) we obtain 



&lon = de 


(3) 


A similar result was arrived at by Chichka 11 using a second-order approximation. 

Equations (1) and (8) show the relationship between the dimension of the relative 
motion in terms of e and i. However, we need to develop an algorithm to place n 
spacecraft in an ellipse-like formation defined by S ton and 6 lat . 
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Figure 3: The Rotating Formation Projected into the ir — x/r Plane 

SPACECRAFT PLACEMENT IN THE ROTATING FORMATION 

In this section an algorithm to place n spacecraft along the ellipse-like path in 
the relative reference frame is developed in terms of 5^ and Si at . The spacecraft can 
be placed in the formation in three different configurations: equal angle separations, 
equal arc length separations, and equal time separations. In this paper an algorithm 
is developed to place the spacecraft around the ellipse-like path with equal time 
separations. For circular formations, equal time separation is equivalent to equal 
arc length and equal angle separation. For ellipse-like formations the relationship 
between the three configurations is non-trivial. The formation is developed by first 
investigating the requirements to place a single spacecraft on an ellipse-like path 
defined by <5^ and 5i at . Then, the method is generalized to include n spacecraft 
spaced equally in time. 

We know from Eqs. (1) and (8) that the dimensions of the relative motion path 
in J- re i depends only on e and i. Therefore, for all the spacecraft to move along an 
ellipse of the same dimensions in J- re h a H the orbits in the formation must have a 
common e and i. Solving for i in Eq. (1) yields 

i = di at /2 (9) 

From Eq. (8) the eccentricity can be written as 

e = Si on /4 ( 10 ) 

To ensure a cohesive formation, all the orbits must also have the same semimajor 
axis. Therefore, only the argument of periapsis. a/, the longitude of ascending node, 
Q, and the true anomaly, u may vary for the orbits in the rotating formation. 
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Figure 4: Time Evolution of the Rotating Formation 


In Figure 4, the position of an arbitrary satellite at two different epochs is shown. 
At the initial epoch the position is described by r fo ; at a later epoch the position is 
described by r tt . The origin of the reference orbital frame is chosen to lie along the 
Xfl-axis at the initial epoch in a circular orbit whose plane of motion is the xr — Ur 
plane. We begin constructing the formation by placing a spacecraft in position 1 
shown in Figure 1. The relative motion is chosen to be clockwise. For clockwise 
motion to occur the spacecraft at position 1 must be moving in the y re i direction 
faster than the origin of the reference orbit. Therefore, spacecraft 1 is positioned at 
the periapsis of its orbit: 

*d = 0 (11) 

The angles o;i and Qi must be chosen to position the periapsis of orbit 1 over i/*. By 
inspection of Figure 4 we see that this is accomplished by choosing 

ill = 37 t / 2 U\ = 7t/2 (12) 


In summary, the orbital elements to position a spacecraft at position 1 in T re .i at 
the initial epoch are 


a[ — a e\ — Si on / 4 — Si at /2 uq — tt/2 Qi — 3tt/2 v i — 0 (13) 


We now wish to generalize the algorithm to place n satellites spaced equally in 
time along the ellipse-tike path seen in Figure 1. For all satellites to move along the 
same path in Fr, all the orbits must share a common a, e, and i. However u,\ Q, and 
v must be determined for the remaining orbits. 
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Tin 1 r i me of separation, t. s , between spacecraft in an equal time formation is given 
by 

t, = T/n = — ./? ( 14 ) 

n y /.i 

To determine the orbital elements of the i th spacecraft in formation orbit 1 is pro- 
pogated through time £, where 

U = (* - l)i, (15) 

The configuration of spacecraft 1 and the reference orbit at the initial epoch, t 0 . and 
at £, can be seen in Figure 4. We need to determine u t> , the value of u for orbit 1 
after the passage of time t t . For small e a first order expansion for u is 

u ti = M ti + 2esin A/ t , (16) 

where 

M ti = M 0 + (i — l)r)t s (17) 

Recall that the true anomaly for orbit 1 at the initial epoch is zero. Therefore, the 
mean anomaly of orbit 1 at the initial epoch, M 0 , is also zero. Substituting Eq. (14) 
and (17) into Eq. (16) the relation for the true anomaly is obtained: 

v ti = — (* — 1) + 2esin ~ 1)^ (18) 


For a Keplerian orbit only v changes with time. Therefore, the orbital elements 
to space n spacecraft equally in time in Tr are known. However, we are interested 
in the elements that yield n spacecraft spaced equally in time in T re i. Therefore, we 
need to determine the new configuration of T re i after the passage of time f t . Then 
the orbital elements to space n spacecraft equally in time in T Te i can be determined 
by a coordinate transformation. 


The frame J- Te i is defined by a circular, equatorial reference orbit in Tr. Therefore, 
the motion of the coordinate system T re i after time ti is simply a 3-rotation through 
the angle 6 = T]t t . Performing a 3-rotation through tjU to T re i is equivalent to rotating 
Q, through an angle 0 = —r/t t . Therefore, to have n spacecraft spaced equally in time 
in T re i the node is written 


o o/o t 37r 

= 3tt/2 - rjt i = — 


(19) 


In summary, to obtain the orbital elements to place the i th spacecraft in an equal 
time, rotating formation in T re i orbit 1 is propogated through f t as defined in Eq. (15). 
Then the longitude of the ascending node is rotated backwards through r}t t to Q t as 
defined by Eq. (19) to account for the motion of Trei- Hence a. e. i. and jj are the 
same for all orbits in the formation. The true anomaly and the longitude of ascending 
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node arc given by Eqs. (18) and (10). In general, the i"' set of orbital elements is 
given by 


", 


" ", = 4,„/o 

Vi 


>■< = tut/- = ?r/2 fi, = 3rr/2 

?(i-l)+2e i sin(a (i _i)) 


J T I - L) 

n 


( 20 ) 


where 1 < i < n. 


PERFORMANCE OF THE ROTATING FORMATION 

We wish to investigate the performance of the rotating formation using an appro- 
priate metric. We propose that the performance of the formation can be evaluated 
by investigating the angular separations between spacecraft over the course of an or- 
bit. 12 Specifically, an orbit-averaged quantity based on the sum of the instantaneous 
separation angles is considered: 

= ~L rf ± «,.(<*) m ( 21 ) 

n > T J o 1 

where w a is an instantaneous position metric to be chosen by the analyst, a t] is the 
angular separation between the i th and j th spacecraft, and tj s is the number of unique 
separation angles given by 

9 

n — n 


A measure of this form allows the analyst to define an instantaneous metric for the 
angular separation between a pair of spacecraft. The performance of all pairs in for- 
mation is evaluated accordingly and summed together to obtain a measure of the 
formation performance at a particular instance in time. By integrating the instanta- 
neous metric over an orbit a measure of orbit effectiveness is obtained. A parabolic 
form is chosen for the instantaneous weight function, w a . In doing so, the flexibility 
is enabled to choose an upper and lower limit for the acceptable angular separations 
between two spacecraft as well as the optimum configuration for two spacecraft. 

To solve the performance problem, the position information for all spacecraft in 
formation over an orbit period is needed. In general, spacecraft are subject to forces 
from the spherical primary as well as perturbations from atmospheric drag, solar radi- 
ation pressure, third body forces, and non-spherical effects. However, approximating 
the orbital dynamics with Keplerian motion is suitable for a preliminary analysis. 
Therefore, to obtain the conditions that yield an optimum configuration Eq. (21) 
must be solved for an optimum where the dynamics of the formation are described 
by the two body equation of motion with the initial conditions given by the orbital 
elements in Eq. (20). 

We now investigate optimum configurations of the formation according to the 
angular separation performance metric. 
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An Analytic Solution 

Wo wish to determine optimal configurations of the rotating formation in terms 
of the orbital elements defining the rotating formation. Hence, we seek the orbital 
elements that yield an optimum in Eq. (21). In general. U Q does not permit a closed- 
form analytic solution. However, some simplifying assumptions are made that allow 
the application of the metric to the rotating formation to obtain analytic results. In 
summary, a closed-form solution for the maximum in Eq. (21) is desired where 

w a (oti j) = CiOtfj + c 2 a tJ + c 3 (23) 


The performance measure seen in Eq. (21) is based on the angular separations 
between spacecraft, which are determined by the orbital elements in Eq. (20). For 
the rotating formation, the orbital elements are simple functions of a, n, 5 lon , and 5 lat . 
However, the angular separation between spacecraft does not change with semimajor 
axis. Therefore, the functional expression for is 

otij = /(nJioruhat) (24) 

The parabolic instantaneous weight function, w a , can be uniquely defined by three 
constants Cj, c 2 , and c 3 . The result is that the performance problem has six design 
variables such that 

W<* = f $loni filati t-l) C 3 ) (25) 

The dimensions 5^ and 5[ at can be written in terms of inclination and eccentricity 
therefore, otij can be expressed in terms of n, e, and i. A solution for a t] in terms of 
n, e, and i allows a solution for W a of the form 

W a = /(n,e,z,ci,c 2 ,c 3 ) (26) 


An expression for can be obtained using the definition of the cross product: 


Tj X Vj ||= TiTj Sin Ctij 


(27) 


For close formations is small so we can assume 


sin otij ~ a i3 


(28) 


Solving for a l} we obtain the expression 


on 


r t x r. 


r x rj 


(29) 


The vector r, in the perifocal system can be expressed as 


a(l-e 2 ) 

r ‘p = 7T 

\ + e cos u, 


cos u L sin i/ t 0 


(30) 
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whore u L and u } an' given hv E(i. (20). The vector r ; can be written using a similar 
expression by replacing the subscript t in Eq. (30) with j. 

To use Eq. (29) the vectors must be rotated into a common reference frame using 
the appropriate rotation matrix. The rotation matrix to go from the perifocal system 
to T r is constructed for the i th position vector using 

Rf p = R 3 (-a)Ri(-0R3(-w) (31) 


where Q, is given in the orbital element set in Eq. (20). Then the inertial represen- 
tations of r t and r j can be written using 


r 


»/ 


R/’r t 


L J/ 


= R /p 




(32) 


Using the preceding, we have the information required to evaluate ot x] from Eq. (29). 
The expression, which is not shown here, is a non-trivial transcendental function of 
the form 

otij = f{n,e,i) (33) 

We now make some assumptions to arrive at a simpler form for W a . 


Simplifying Assumptions 

In general, solving W a requires integrating the motion over an orbit. However, 
for special configurations of the rotating formation, the relative spacing between the 
spacecraft is nearly constant over an orbit. This is true for near circular rotating 
formations when 

~ 5i at (34) 

The relative positions are also nearly constant over an orbit for formations com- 
posed of a large number of spacecraft. In Figure 5, a rotating formation of four 
spacecraft is shown. The shaded circles represent the positions of four spacecraft at 
the initial epoch. The white circles represent the positions of the spacecraft after 
a small passage of time, At. The relative positions, and hence the angular separa- 
tions, have changed after the time At. However, when the spacecraft at position 1 
reaches position 2 the relative positions of the spacecraft will be nearly the same as 
at the initial epoch. The relative positions are not exactly the same because all of the 
spacecraft in formation are at different points in elliptical orbits. For example, for 
this configuration spacecraft 2 is at periapsis and spacecraft 1 is at a true anomaly 
given by Eq. (20) with i = 4. Therefore, the time for spacecraft 1 to reach position 2 
is slightly different from the amount of time for spacecraft 2 to reach position 3. For a 
four-spacecraft formation the relative positions are nearly the same after a passage of 
time T /4. In general, for a formation of n spacecraft, the relative positions are nearly 
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Figure 5: Rotating Formation of Four Spacecraft 



Time, s 


Figure 6: Integrand vs. Time for a Highly Elliptical Formation 


equivalent after time T jn. Therefore, the period of the formation decreases with n. 
More importantly though, as n increases the amplitude of the integrand in Eq. (21) 
goes to zero. This is demonstrated in Figure 6. A plot of the integrand over an orbit 
period for n = 4, 6, 10 is shown for a highly elliptical formation with d^ = 6000 m 
and di at = 1000 m. As n increases the period and amplitude of the motion decrease. 
Hence, the relative positions of the spacecraft are nearly constant for large n. 


If the relative motion of a formation is nearly constant over an orbit the integral 
can be neglected in 

W a = 4 ; / E E w(oij) dt (35) 

VsT Jo fr{ - =l+1 
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and U** can bo expressed as \\' n where 


n - 1 n 



(36) 


With Eqs. (23), (29), and (36) the tools are available to express \V a analytically. 
We are interested in the dimensions of the rotating formation that result in optimum 
performance. Recall that the shape of the rotating formation is determined by e and 
i. The analytic expression of \V a is a non-trivial function in terms of six variables. 
A simplification of the expression is necessary to solve for optimum configurations in 
terms of e and i. Recall that W a is not dependent on a. Therefore, the analysis can 
be continued without concern for the specific value of semimajor axis. 

To obtain an expression for W a , we need to evaluate the summation in Eq. (36) 
for a given n. A value of n = 4 is chosen to demonstrate the techniques used to obtain 
the optimum. Results for several values of n are presented later to determine trends 
in the optimum solutions for changing n. 

Assuming n = 4 in Eqs. (23), (29), and (36), W a can be expressed in terms of e, i, 
Ci, C 2 , and C 3 . The results can be expanded to third-order in small e and i about e = 0 
and i = 0 to obtain an approximate expression. However, the resulting expression 
has a singularity at e = 0. Plotting the exact expression for W a in negative i and e, 
as shown in Figure 7, lends some understanding: there is a cusp at e = 0 and i — 0. 



Figure 7: The Singularity in Eccentricity 
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Tin* singularity in r can bt' avoided by expanding about, another point. To do so 
we make a substitution using 

e = ^(1 + e) (37) 

where 0 < f < 1 is a small dimensionless parameter. Expanding the resulting 
expression to third-order in i and e yields an approximate expression for Id a : 

... / 1 lc 2 13c 2 \ i Scii 2 e 8 c { i 2 c 2 e 2 i ( n/2c 2 4- c 2 

W - = - (lF + •' + — ■ + — ■ 3 

(38) 

To solve for the e and i which results in an optimum solution, we must solve the set 
of simultaneous equations 


\ 2n/2c 3 4- 2c 2 

j«+ j +c, 


aw. 

di 


= 0 



(39) 


for e and z. However, expressions for dWJdi and dW a /de in terms of e must be 
found. Differentiating Eq. (38) with respect to i yields 


dW a _ 64d(l 4- e)i + c 2 (8(l + y/2) + 4(1 + y/2)e + v/2e 2 - (22 + 13y/2)z 2 
di ~ 12 


To obtain an expression in terms of e a substitution using 



i 


(40) 


(41) 


is performed to obtain 

dW a _ 128 Cl ez 2 4- c 2 (4y/2e 2 + 4(2 + y/2)ei + i 2 ( 4 4- 5%/2 - (22 + 13y/2 )i 2 )) , 

di ~ 12z 2 ( } 

Taking the partial derivative of W a with respect to e yields 

dW a i 2 (c 2(2 T 2y/2 4- \/2e) + 16ciz) ^ ^ 

nr = 12 (43) 

To determine an expression in terms of e a multiplication by de/de is performed where 


The derivative with respect to e is 


de_ 

de 


= 2 fi 


dW a _ dW a de _ * (l6c L z 2 + c 2 (2s/2e + (2 + 
de de de 12 


(44) 


(45) 


Solving Eqs. (42) and (45) simultaneously yields expressions for e and i that yield an 
optimum in W a . The expressions are complex and not shown here. However, there 
is a relation between the optimum e and i such that 


^ ’opt 


- i opt 1 2 


(46) 
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Recall that we assume n = 4 in order to obtain an explicit expression tor ll',. 
Therefore, E<|. (46) appears to be dependent, on the assumption that the forma- 
tion is composed of four spacecraft. However, the same result is obtained using 

n = 2, 3, 4. 5. 6. 8, 10. 


Equation (46) has several interesting physical implications. Substituting the re- 
lations for e and i from Eq. (20) into Eq. (46) yields 6i on = 8i at . Hence, the ideal 
rotating formation according to lV a is circular. 


By substituting e = i/2 into Eq. (36) an expression is obtained for U Q in terms of 
n, i, Ci, c 2 , and C 3 . Recall that the only assumptions made to arrive at Eq. (36) are 
small angles and that the integrand in Eq. (35) is constant in time. To simplify W a , 
we expand in a power series to second-order in i and investigate a larger formation 
with n = 8 . A second order expansion in i for n = 8 yields 


W a = ^ 


4 2^2 2 

c 2 I 7 = 3 = b - | 1 + C 3 

Jyj 2-y/2 7 7 


(47) 


To provide a more intuitive understanding of the results, a new form for the parabolic 
weight function of Eq. (23) is introduced: 


w Q ( c%ij ) — (ouj ) {cx t j c?i ) (48) 

where a u and on are the upper and lower limits on angular separation respectively 
and Ofj is a normalization factor. The constants C\, c 2 , and c 3 from Eq. (23) can be 
written in terms of a u , a t , and a s using 


Cl 

= 1 /<*s 

(49) 

C 2 

= -(a u + ai)/a s 

(50) 

C 3 

S' 

3 

0 

II 

(51) 


For circular formations, the radius of the formation is simply the inclination of 
the orbits composing the formation. Therefore, solving for the optimum radius, r 0 , 
from Eq. (47) and substituting for ci, c 2 , and c 3 using Eqs. (49-51) we obtain 

r o = fjg + 4- — {\fz — \/2 4- \fl 4- n/2^ (an 4- a u ) (52) 

A new angle, a m , is defined such that 


& m — 




(53) 


This angle is the ideal separation angle according to the instantaneous weight function 
w a . The optimum radius, r a , can be rewritten in terms of a m : 


r„ = 


'1 

8 + 






(54) 
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The samt* techniques ran he appled to find tin* optimum radius fur different values 
of n. Tin* results for different values of a are 
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10 

71 = 12 

r o 



1 




12 


'10 v/5 /l0 ^5 


Otr. 


— — ^2 + V2 + \/3 + (^2 — \/3 + \fl + V3) ''j oc r , 


( 00) 

(56) 

(57) 

(58) 

(59) 

(60) 

(61) 

(62) 


There are several interesting trends seen in the equations for r Q using different values 
of n. The equations are functions only of a m . Therefore, the optimum configuration 
for a given n depends only on the ideal separation angle according to the instanta- 
neous weight function w a . The shape of the parabola does not affect the results. 
Furthermore, with the exception of n = 3, the optimum radii solutions are linear 
functions of a m . 


The relation between r a and n is not immediately obvious. However, upon close 
inspection we can determine a more general form of r D for several values of n. For 
n — 2,4,6 r 0 can be expressed as 

r ° = n (~2~ + ^>) am (63) 

A plot of the optimum radius over a range of n values is seen in Figure 8. The 
points are the analytic solutions according to Eqs. (55-62) for a m = 0.000375 rad. 
The curve is a plot of the empirical solution in Eq. (63). The points for n = 2,4,6 
fall exactly on the curve. This result is expected because the empirical function is 
derived from the equations of r a for n = 2, 4, 6. The values of r a for n = 3,5 are in 
good agreement with the empirical solution. However, for n > 8 the accuracy of the 
empirical function degrades. The empirical function predicts a downward trend in r a 
for n > 8. However, the analytic solution appears to approach a limiting value. 

The limit of r 0 as n approaches infinity can be determined using simple geometric 
relations. For n — >• oc, w p (a i; ) is the same for all /. where 1 ^ j. Therefore, we only 
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Figure 8: Optimum Radius for the Rotating Formation 
need to evaluate w p {ctij) for an arbitrary i. Thus for n — ► oo 

= ri-^+oo V E E = J™ ~1 E W p( a j) (64) 

'» i=l j=n+l 1 1 j=2 

Figure 9 shows two arbitrary radius vectors, rq and r 2 , defining two spacecraft posi- 
tions in a circular formation. Angle 77 is the angular radius of the formation. The 
angle around the circle between iq and r 2 is denoted 9. Using these definitions Eq. (64) 
can be written as 


1 n If 2 * 

W *> = Jim E w p( a i) = 2^ L w p^ dd ( 65 ) 

j = 2 

From inspection of Figure 9 we can write 

r x = [sin 77 0 cos rf] T (66) 

r 2 = [sin 77 cos 9 sinry-sin# cosr/] T (67) 

Recall that W a is not dependent on the semimajor axis. Therefore, we use unit 
vector representations of the spacecraft positions without loss of generality. Using 
ri = r 2 — l, a can be expressed as 


a =11 ri x r 2 || 

for small a. Substituting Eqs. (66-67) into Eq. (68) we obtain 

a = Zsin 2 77 ( — 2 cos 2 r f { — 1 + cos 9) + sin 2 77 sin 2 9^ 


( 68 ) 


(69) 
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Figure 9: Circular Formation Geometry for n = oo 


Substituting Eqs. (69) and (53) into Eq. (48) the integral in Eq. (65) can be evaluated: 

1 


= 


167ra, 


(—7 7T — 167T ona u + 37rcos4r^ + 8a m (A;2 — A^)) 


+ (4cos2r / (tt + 2a m (k 2 - h)) + Z2a m \/sin* 77 ) 

where k\ and k 2 are given by 


(70) 


k i 

k 2 


In 

In 


i-/ 


sin 2 t 7 


1 4- ysi 


sin 2 r f 


( 71 ) 

(72) 


To solve for the optimum radius as n —¥ oo we need to solve 


dWoo 

drj 


= 0 


(73) 


However, the resulting derivative is a transcendental function in terms of a m and 
Tj and cannot be explicitly solved for the optimum radius. Therefore, Eq. (70) is 
expanded in small r f to second order: 


W n 


1 ( 2 rj _ 


a 


7 T 


+ Oil Ot 




(74) 


Defining as the optimum radius as n — ► oc and solving Eq. (74) for obtains 



7T 


(75) 
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* 2 nd Order Solution 

Function Fit 

Asymptote 

a m -0.000375 

n*4 


2 4 6 8 10 12 

Number of Satellites in Formation 

Figure 10: Optimum Radius with Asymptote vs. n 

Like r 0 for small n, is only a function of a m . Therefore, the shape of the 
parabolic instantaneous weight function does not affect the optimum solution for 
large n. In Figure 10, the solution of r a for a range of n values is plotted with the 
asymptote r^. The approximation for the asymptote is in good agreement with the 
analytic solution for r 0 for large n. 

Recall that several simplifying assumptions are made to obtain an analytic solution 
for W Q and the limiting value as n — > oo. We need to determine if the simplifications 
yield solutions that are in good agreement with the exact solutions. Plots of W a for 
n = 3, 4, 8, and 12 are shown in Figures 11-12 for a range of a m . The data points are 
the exact numerical solutions for different radii obtained using a Runge-Kutta fourth 
order integrator to solve Eq.(21). The curves are the approximate analytic solutions 
for W a . The approximate solutions are in excellent agreement with the exact solution 
for low a m . For high a m , on the order of a m = 0.25, accuracy of the approximate 
solution for W a is poor. A decrease in accuracy is expected for high a m because the 
small angle approximation was used in the development of the approximate solution. 
However, the optimum in performance occurs at nearly the same point according to 
both the exact and approximate solutions even for ot m = 0.25. In Figure 13, plots 
of the exact vs. approximate optima are shown for a range of a m . The approximate 
solutions are accurate even for large a m on the order of a m = 0.5. In the range 
of 0.5 < a m < 1 the accuracy of the approximate solution degrades. However, the 
approximate results are still useful in this range. We suggest that for 0.5 < a m < 1 
the approximate solutions be used as an initial guess to solve for the exact optimum 
using numerical techniques. 




Optimum Radius, rad Waighi Function Value W«^N Fui^bon Vak* 





Figure 13: Approximate and Exact Optima vs. a m 
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CONCLUSIONS 


In this paper an algorithm is developed to place n satellites equally spaced in 
time around an ellipse in the relative reference frame in terms of the classical orbital 
elements. The dimensions of the ellipse-like path are shown to be simple functions of 
orbit eccentricity and inclination. The performance of the formation is investigated 
using a performance metric based on the orbit averaged sum of the angular sepa- 
rations between spacecraft. The exact expression is a non-trivial function of orbit 
eccentricity, inclination, the number of spacecraft, and the constants that define the 
instantaneous weight function. However, several simplifying assumptions allow for an 
analytic approximation. Using the resulting approximate equation for VU Q , optimal 
configurations are found for different values of n. We show that, according to the angu- 
lar performance measure, the optimal rotating formation is circular. The optimal radii 
for different values of n are found to be simple functions of a m , the ideal separation 
angle according to the instantaneous weight function. Interestingly, the optimal con- 
figuration does not depend on the shape of the parabola chosen for the ins tant aneous 
separation metric. For n < 8, the optimal radius is r a = ((n — 2)/2 -I- \Jn/2)a Tn /n, 
whereas for large n (> 8), the asymptotic limit r 0 — ► 2a m /7r is quite accurate. 


ACKNOWLEDGMENTS 

The first author was supported by a grant from the Guidance Navigation and Control 
Center at NASA Goddard Space Flight Center, under the direction of David Folta 
and David Wiedow. The second author was supported by a grant from the Air Force 
Office of Scientific Research, under the direction of Arje Nachman. 


NOTATION 

a semimajor axis, m 

Ci,C 2 ,C 3 constants defining parabolic weight 
d spatial separation between two spacecraft, m 

e orbit eccentricity 

Ti inertial reference frame 

JF re / relative reference frame 

Tr rotating reference frame 

Tp perifocal reference frame 

i orbit inclination, rad 

l true longitude at epoch, rad 

n number of spacecraft in formation 

n s number of distinct angles in a formation 

r x optimum radius as n — > oo . rad 

r„ optimum radius for rotating formation, rad 
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T orbital period, s 

t, propagation time for i' h spacecraft, s 

t. s e([ual time separation for n spacecraft, s 

orbit performance measure 
H x U" rt as n — > oc 

VV a approximate expression for lV a 

w a instantaneous position metric 

a t] angle between the i th and j th satellites, rad 
at lower limit on angular separation, rad 

ct m ideal angular separation, rad 

a s scaling factor, rad 

a u upper limit on angular separation, rad 

Sion formation longitudinal separation, rad 

Si at formation latitudinal separation, rad 

77 orbit mean motion, rad/s 

v true anomaly, rad 

Q longitude of ascending node, rad 

u argument of periapsis, rad 
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